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On Log-concavity of the Generalized Marcum 
Q Function 
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Abstract — It is shown that, if v > 1/2 then the generaUzed 
Marcum Q function Q^{a,b) is log-concave in & G [0, oo). This 
proves a conjecture of Sun, Baricz and Zhou (2010). We also 
point out relevant results in the statistics literature. 

Index Terms — increasing failure rate; log-concavity; modified 
Bessel function; noncentral chi square. 



I. Introduction 

The generalized Marcum Q function 1141 has important ap- 
plications in radar detection and communications over fading 
channels and has received much attention; see, e.g., |[3l, lIU, 
ifTOl . |[T3]-|[17J and li9J-l21J. It is defined as 



Quia,b) 



T exp 



I^-i{at)dt (1) 



where > 0, a,b > and denotes the modified Bessel 
function of the first kind of order i' defined by the series |jjj 
(9.6.10) 

(t/2)2fc+'' 



k=0 



kir{iy + k + l)' 



(Q^{0,b) is defined by taking a I 0.) Recently, Sun, Baricz 
and Zhou flT] have studied the monotonicity, log-concavity, 
and tight bounds of Q^{a, b) in great detail. We are concerned 
with log-concavity, which has intrinsic interest, and can help 
establish useful bounds; see |21 ) and the references therein for 
the large literature in information theory and communications 
on numerical calculations of Q^{a,b). 

This note resolves some of the conjectures made by f2T\. 
We also point out relevant literature in statistics on both 
theoretical properties and numerical computation of Qiy{a, b). 
Our Theorem [U proves Conjecture 1 of [21]. 

Theorem 1: The function Q,y{a,b) is log-concave in b E 
[0, oo) for all a > if and only if > 1/2. 

A sufficient condition for log-concavity of an integral like 
(HJ is that the integrand is log-concave in t. Proposition [T] and 
Theorem |2] take this approach. 

Proposition 1: The integrand in ([T]l is log-concave in t E 
(0, oo) for all > 1/2 if and only if < a < 1. 

Theorem 2: The integrand in ([TJ is log-concave in t E 
(0, oo) for all a > if and only if v > vq where vq k, 
0.78449776 is the unique solution of the equation 



/,(V5^^) 



2v 



in the interval v E (1/2,3/2). 

Note the difference between Proposition [T] and Theorem |2] 
the former gives a criterion for log-concavity in t for all 
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J/ > 1/2 whereas the latter gives one for all a > 0. From 
Proposition [T| and Theorem |2] we obtain Corollary [T] which 
confirms part of Conjecture 2 of [21]. 

Corollary 1: The function 1 — Q^(a,b) is log-concave in 
b E [0, oo), if either (i) > 1/2 and < a < 1, or (ii) > i^o 
as in Theorem |2] 

The case of Qi{a,b) (Marcum's original Q function) is 
especially interesting. If = 1 then the integrand in ([TJ is 
the probability density function (PDF) of a Rice distribution, 
Qi{a,b) being the corresponding tail probability, or survival 
function. Therefore Theorem |2] yields 

Corollary 2: The probability density function, cumulative 
distribution function (CDF), and survival function of a Rice 
distribution are all log-concave. 

In general, let X be a noncentral random variable with 
21/ degrees of freedom and noncentrality parameter a^. Then 



g^(a, 6) =Pr (VX > 6 



Equivalently, 1 — Q^{^/a, \/b) is the CDF of a noncentral 
random variable with 2v degrees of freedom and noncentrality 
parameter a. The noncentral distribution plays an important 
role in statistical hypothesis testing and has been extensively 
studied. We mention ||6l, lfT2l on numerical computation and 
Q, in, lITSl on theoretical properties. Its CDF, and hence 
Qv{a,b), can be routinely calculated (e.g., using pchisq() in 
the R package). 

Concerning theoretical properties, Finner and Roters [7^ (see 
also [5J) have obtained the following results using tools from 
total positivity ifTH . 

Theorem 3 (fl), Theorems 3.4, 3.9; Remark 3.6): The 
function 1 — Q^{^/a, \/b) is log-concave 

. in 6 e [0, oo) for > 0, a > 0; 

• in v > Q for a, b > 0; 

. in a > for > 0, 6 > 0. 
The function Qu{\/a, Vb) is log-concave 

• in 6 e [0, oo) for > I, a > 0; 
» in ly E [1/2, oo) for a, b > 0; 

. in a > for > 0, 6 > 0. 
Theorem [3] and Corollary [T] cover several results of lETI . 
including part of their Conjectures 2 and 3 (see also ||20) ). 
The parts of these conjectures that remain open are 

• l—(5;y(a, 6) is log-concave in 6 e [O,oo)fori^e [1/2,1/0) 
and a > 1; 

• Qiy{a, b) is log-concave in 1/ e (0, 1/2] for a, b >Q. 

In Section II we prove Theorems [T] |2] and Proposition [T] 
The proof of Theorem [T] uses a general technique which may 
be helpful in related problems. The proof of Theorem |2] relies 
partly on numerical verification as theoretical analysis appears 
quite cumbersome. 

II. Proof of Main Results 

The following observation, which is of independent interest, 
is key to our proof of Theorem [1] 

Lemma 1: Let f{t) be a probability density function on 
R = (—00,00). Assume (i) f{t) is unimodal, i.e., there exists 
to E H such that f{t) increases on (—00, to] and decreases 



2 



on [to,oo); (ii) /(to-) < f{to+); (iii) f{t) is log-concave in 
the declining phase t e {to,oo). Then the survival function 
F{b) = Jj^ f{t) dt is log-concave in & e R. 

Proof: Assumption (iii) implies that F{b) is log-concave 
in & € [io,oo). Because f{t) increases on {—oo,to], we 
know F{b) is concave and hence log-concave on (— oo,to]- 
By Assumption (ii) we have 

F'ito-) = -f{to-) > -f{to+) = F'{to+). 

Hence F{b) is log-concave in 6 e R overall. ■ 
Remark 1. A distribution whose survival function is log- 
concave is said to have an increasing failure rate (IFR) [4]. 
Distributions with IFR form an important class in reliability 
and survival analysis. Lemma [T] provides a simple sufficient 
condition for IFR distributions. 

Henceforth let /(<) be the integrand in ([T]l for t > 0. 
Equivalently, f{t) is the density function of a noncentral x 
random variable with 2^ degrees of freedom. Define 



I At) 



(2) 



We use r'jj{t) to denote the derivative with respect to t. 

Lemma 2: If > 1/2 then f'{t)/{tf{t)) decreases in t e 
(0,(X3). 

Proof: Let us assume v > 1/2 and a > 0. The boundary 
cases follow by taking limits. Direct calculation yields 

nt) _ y ^ aK-M) 
tf{t) tl,-i{at) 



2v — 1 ar^{at) 

where (|3]l uses (|2]l and the formula [ 1 1 (9.6.26) 

/:_i(i)=/.(<) + ^/._iW. 



(3) 



(4) 



Since (2i^ — decreases in t, we only need to show that 

Ti, (t) 1 1 decreases in t. We may use the integral formula of [ T\ 
(9.6.18) and obtain 

t i2,.-l)J^gis,t)ds 
g{s,t) = {1- s^y-^^^coshits). 



where 



As can be easily verified, if < ti < t2 then g{s, t2)/g{s, ti) 
increases in s G (0,1). That is, g{s,t) is TP2 HI]. Since 
1 — decreases in s e (0, 1), by Proposition 3.1 in Chapter 1 
of ifTTI . the ratio (1 — 5^)3(5, t) ds/ g{s,t) ds decreases 
in i G (0,00), as required. ■ 
Proof of TheoremU] Let us assume 1/ > 1/2 and show 
log-concavity. By Lemma |2l either (i) f'{t) < for all t E 
(0, 00) or (ii) there exists some to e (0, 00) such that f'{t) > 
when t <to and f'{t) < when t > to- (Since f{t) dt = 
Qi/(a,0) — 1, it cannot happen that f'{t) > for all t E 
(0,00).) In either case /(<) satisfies Assumptions (i) and (ii) 
of Lemma [T] (/(i) = for t < 0). Let us consider Case (ii); 



the same argument applies to Case (i). For t e (to, 00) we 
have f'{t) < 0, and hence 

1 d^ 1 f'(t) 



log fit)] <0 



1 d 

dt yldt 

where the last step holds by Lemma |2] Thus f{t) is log- 
concave in t G (<o,oo) and Assumption (iii) of Lemma [T| 
is satisfied. We conclude that Qi^{a,b) = f{t)dt is log- 
concave in b G [0,00). 

It remains to show that, if Q^{a,b) is log-concave in 6 e 
[0,00) for all a > 0, then we must have ly > 1/2. Let us 
consider a ~ 0. We have 

Q.(0,6)^l-^^%-e-/^d. 

As 6 I 0, it is easy to see that log Qi,{0, b) behaves like 

log (1 - Cb^" + 0(52-)) = -Cb^'^ + o{b^'') 

with C = 2-''/T{v + 1). Hence, if < 1/2 then Q^(0, b) is 
no longer log-concave for b near zero. It follows that the 1/2 
in Theorem [T] is the best possible. ■ 
Proof of PropositionU} Using (|3j we get 



^logfit) ^ -1 + aXiat). 



(5) 



However, 



= 1 —r.it)~rlit) 



(6) 



where (|6]l holds by applying (|2]i, (HJi and the recursion ([T] 
(9.6.26) 

i.+i{t) = i.-i{t)~ ^Ut). 

If > 1/2 and < a < 1 then rl{at) < 1 hy and we 
have 

^log/(t)<a2-l<0. 

Hence f{t) is log-concave in t e (0,oo). 

To show the converse, suppose f{t) is log-concave in t for 
all 1/ > 1/2. Consider = 1/2. As t J, we have r^{t) 0, 
and d^ log f{t)/dt^ — ^ — 1. Hence we must have a < 1. ■ 

Remark 2. For 1/ > 1/2, the function f{t) is log-concave 
in its declining phase, as shown in the proof of Theorem [T] If 
a G [0, 1] in addition, then Proposition[T]shows that f{t) is log- 
concave in all t e (0,c»). For a > 1 and v > 1/2, however, 
numerical evidence suggests that f{t) may not be log-concave 
in its rising phase. Hence a version of Lemma [T] cannot be 
applied to 1 — Q^{a,b). Log-concavity of 1 — Qiy{a,b) in b 
appears to be a difficult problem. 

Let us establish two lemmas before proving Theorem |2] 

Lemma 3: The function f{t) is log-concave in t G (0, 00) 
for all a > if and only if the function 

Kit) = 1 ^ - ^r.(t) - rlit) (7) 



f2 



t 
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is nonpositive for t e (0, oo). 
Proof: By (|6]l we get 

Kit) = r'Jt) - 

If h^{t) < then by © we have 

d2 



(8) 



dt2 



log f{t) = a^K{at) - 1 < 0. 



Conversely, if f{t) is log-concave in t G (0, oo) for all a > 0, 
then holding at constant while letting o — t- oo yields /i^(s) < 
for each s € (0, oo). ■ 
Lemma 4: The function 

strictly increases in v E [1/2,3/2] and has a zero at i^o « 
0.78449776. 

Proof: Although this only involves a one-variable func- 
tion over a small interval, it is verified by numerical calcula- 
tions, as theoretical analysis becomes complicated. The value 
of i/Q is computed by a fixed point algorithm. ■ 

Proof of Theorem |2} Define hi, (t) as in (|7|i and h'o as 
in Lemma|4] We examine the intervals (0, 1/2], (1/2, i^o) and 
[I'OjOo) for ly in turn. If0<i^<l/2 then letting i | we 
have r^{t) -> and h^{t) > for small t. By Lemma [3] f{t) 
is not log-concave for all a > 0. 

Let us assume > 1/2. Differentiating (|7]l with respect to 
t and applying ([8]l we get 



Kit) = - 



2jy - 1 



lu{t) - 



+ 2r,{t)]h,{t) (9) 



where 



L{t)=r,{t)- 



3-2iy 
t ' 



(10) 



For ly > 1/2 we know r,^{t) increases from to 1 as < 
increases from to oo (see 121). Hence, if 1/2 < < 3/2, then 
lv{t) strictly increases and /^(t) = has a unique solution, 
say at ti G (0,oo). If 1/2 < v < i^q, then by Lemma |4] 
lu{V^^-2v) < 0, and hence ti > — 2i/. In view of (Q 
and ( fTOb we have 

, , , 21/ -1 2iy-l/3~2iy\ (3-2iy)^ 



5-21^ 



'-1 



> 0. 



(11) 
(12) 



By Lemma [3] /(t) is no longer log-concave for all a > 0. 

Suppose ly > j/q. We have h,y{t) — > — oo as i ^ and 
h,^{t) — > as i oo. If hi^{t) does become positive, then 
there exists a finite to > such that h,y{to) = and /i^io) ^ 
(at least one sign change should be from — to +). We get 
lAto) < from ^.lfu> 3/2 then ^ yields l^{to) > 
1^1^ (to) > 0, a contradiction. Hence h^{t) < for all t £ 
(0,oo) if 1/ > 3/2. 

Suppose va < V < 3/2. If li,{tQ) — hl{to) = then we 
deduce to = \/5 — 2j/ from O and (|9]l by a calculation similar 
to (fTTTi-lfnii. But /y(V5 - 2iy) = contradicts Lemma |4] 
Hence we may assume ft.'y(to) > and liy{to) < 0. By 



Lemma ID we have ~ 2v) > 0. Because li,{t) is strictly 

increasing, and ti is the solution of li,{t) = 0, we obtain 
to < ti < \/5 — 2iy. The calculation (fTTI) - (fT2l) now yields 
hi^iti) < 0. Because h^{to) = 0, h'^{to) > there exists 
e (io,ii) such that /i^(t*) and /i^t*) < 0. By Q, 
we get Z^(t,) > 0, which contradicts the strict monotonicity 
of l,y{t) as l^{ti) = 0. It follows that h^{t) < 0, t € (0,oo), 
and f{t) is log-concave. Taking the limit we extend this log- 
concavity lo ly = 1^0- ■ 
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